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ABSTRACT 

We consider a recently-proposed alternative explanation of the CV period gap in terms 
of a revised mass-radius relation for the lower main sequence. We show that no such 
thermal-equilibrium relation is likely to produce a true gap. Using population syn- 
thesis techniques we calculate a model population that obeys the claimed equilibrium 
mass-radius relation. A theoretical period histogram obtained from this population 
shows two prominent period spikes rather than a gap. We consider also recent argu- 
ments suggesting that the period gap itself may not be real. We argue that, far from 
demonstrating a weakness of the interrupted-braking picture, the fact that most CV 
subtypes prefer one side of the gap or the other is actually an expected consequence 
of it. 
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1 INTRODUCTION 

The dearth of cataclysmic variables (CVs) between ~ 2 and 
~ 3 hr is well known as the CV period gap. The conven- 
tional explanation of this (see e.g. King, 1988; Kolb, 1996; 
Ritter, 1996 for reviews) invokes the idea that mass trans- 
fer on a timescale shorter than the secondary's thermal time 
will cause this star to be larger than its thermal-equilibrium 
radius. If there is a sudden drop in the orbital braking once 
the secondary becomes fully convective, it is plausible that 
a gap of the observed width results. The lack of a cogent 
predictive theory of the braking mechanism for CVs above 
the gap means that no explanation of this kind can be to- 
tally compelling. However, the interrupted-braking picture 
does naturally explain several features such as the sharpness 
of the gap sides, and (to some extent) the distributions of 
various CV subtypes across it. It has shown remarkable re- 
silience in withstanding various alternative suggestions (e.g. 
Livio & Pringle, 1994 - see King & Kolb, 1995). 

In a recent paper, Clemens et al. (1998; hereafter C98) 
suggest another alternative theory for the formation of the 
gap. This is based on a proposed revision of the mass-radius 
relation for the lower main sequence, as derived from colour- 
magnitude diagrams of nearby stars (Reid & Gizis, 1997). 
C98 derive radii and masses using bolometric corrections 
and the temperature scale from Leggett et al. (1996), and 
an empirical mass-magnitude relation (Henry & McCarthy, 
1993). C98 then assume that CV secondaries do not signif- 
icantly deviate from this equilibrium mass-radius relation. 
Thus the effective mass-radius index 



dln_R2 
din A/a 



(1) 



(-R2, M2 are the donor mass and radius) along the evolu- 
tion above and through the period gap region must always 
equal the local equilibrium index of the main sequence. In 
C98's interpretation, ^ takes values close to 1/3 in two dif- 
ferent sections of the main sequence. Any range of secondary 
masses M2 with ( — 1/3 maps on to a single binary period 
(see eqs. (^, (^ below), in this case to about 3.4 hr and 
2.0 hr. C98 show that this gives a rather slow evolution in 
period at these two points, and claim that this reproduces 
the CV period gap. 

We shall not comment on C98's interpretation of these 
data, except to say that it is certainly not the only possible 
one. Our main point is to take issue with C98's conclusion 
that their revision allows an alternative explanation for the 
period gap. We shall show that, far from producing a period 
gap, this revision in fact produces a pair of period 'spikes' 
at 3.4 hr and 2.0 hr. The CV discovery probability between 
these two spikes is no lower than outside them. We derive 
the condition for a mass-radius relation to produce a gap, 
and point out that this is very unlikely to be satisfied by any 
thermal-equilibrium mass-radius relation. Finally we briefly 
consider the argument by Verbunt (1997) questioning the 
reality of the period gap. 
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2 SECULAR EVOLUTION WITH A FIXED 
MASS-RADIUS RELATION 

Roche geometry allows one to express the orbital period 
change in terms of The secondary's Roche lobe radius 
can be approximated as 



1/3 



(2) 



(Paczyriski, 1971), where M is the total binary mass, C ~ 
0.462 and a is the binary separation. Combining this with 
Kepler's 3rd law shows that the binary period is given by 



From we have 



R2 _ M2 
hence with 



for the rate of change of period 



P 
P 



3C - 1 M2 
2 M2' 



(4) 



(5) 



We thus see that P vanishes where ^ = 1/3. 

The binary evolves in response to the loss of orbital 
angular momentum J. This may occur through gravitational 
radiation or other processes, e.g. magnetic braking. With 
primary mass AIi we can write 
1/2 



J = M1M2 



Ga 
A? 



(6) 



In the following we shall assume that all the mass lost by 
the secondary is accreted by the primary, so that the total 
binary mass remains fixed, i.e. Afi — —M2 > 0, Af — 0. 
Then logarithmically differentiating M n) gives 



Rl 



2J_ 
J 



2M2 ( 5 _ M2 
M2 \&~ Mr 



(7) 



For stationary mass transfer (i.e. M2 — 0) the Roche lobe 
and stellar radius must move in step, i.e. _Rl/^l = R2/R2- 
With (^) and (Q) this gives the mass transfer rate as 



M2 
M2 



j/J 



+ M2/M1' 

The denominator of this equation is automatically positive 
if we require mass transfer stability, i.e. that in the absence 
of angular momentum losses the Roche lobe would expand 
away from the star [Rh/Rh > R2/R2) in response to mass 
transfer. This follows from setting j = in and com- 
paring with and in general requires that the mass ratio 
M2/M1 should be ^ 1. 

Thus, specifying J in some way, the binary evolution 
is fixed. In particular we can express the mass transfer rate 
—M2 and period derivative P as functions of the secondary 
mass M2 ■ The probability p(log P) of discovering the system 
in a given range of log P will vary as 



p(logP) oc 



M2 



\P/P\ 



(9) 



where a is some positive power. For a bolometrically flux- 
limited sample drawn from an isotropically distributed pop- 
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Figure 1. Period histogram p(IogP) (hatched, scale on the left) 
predicted by the proposed mass-radius relation of C98, with the 
cumulative distribution superimposed (solid curve, scale on the 
right), for CVs with P < 5 hr. The figure assumes a = 1 in 
eqn. M), as appropriate for a (visual) magnitude-limited sample, 
and was obtained with —M2 and P/P taken from C98. 



ulation we have a = 3/2. Power-law fits to the bolomet- 
ric correction show that a ~ 1 is appropriate for a visual 
magnitude-limited (mvis-limited) sample (see below), and 
indeed this value gives a better fit to the observed distribu- 
tion in the conventional picture (Kolb 1996). 
From (H, m we get 



p(logP) oc ATa" 



(3C - 1) 



M2 
Ah 



7: + i;-Tf ■ (10) 



For the particular choices a = 3/2 and 1 we have 



p(logP) oc M. 
and 



3/2 



1/2 



(3C- 
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M2 
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p(logP) oc 



M2 
3C-: 



(11) 



(12) 



3 THE THEORETICAL PERIOD HISTOGRAM 

The last three equations show us the effect of any choice 
({M2) on the discovery probability, which is the quantity 
predicting the CV period histogram. Since in general A/2 
and I J/ J I are smooth functions through the evolution, the 
character of the predicted histogram is controlled by ^. The 
choice of C{M2) given by C98 has ^ ~ 1/3 for secondary mass 
ranges which translate to periods 3.4 and 2.0 hr. But from 
( p^ we see that this actually makes p ^ 00 at these two 
points rather than reducing p towards zero between them, as 
required for a period gap. Fig. ^ shows p(log P) for the case 
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a = 1, evaluated directly from with —M2 and P/P from 
C98. As can be seen, this does not give a period gap at all, 
but two spikes. The small period derivative at P = 3.4 hr 
raises the discovery probability dramatically there, as the 
systems spend a long time near this period without being 
particularly faint there. A similar but weaker effect is seen 
at P = 2.0 hr. Note that the discovery probability between 
the two spikes is not lower than that outside them. The case 
a = 3/2 (eq. ^) gives similar results, but with even larger 
spikes. Neither of these distributions looks at all like the 
observed one (see Fig. |^ below), which has roughly equal 
numbers above and below the gap, with approximately uni- 
form distributions in log P in each case. The disagreement 
with observation is even worse for a = 3/2. 

The discovery probability (^) refers to a particular CV 
system at period P with current mass transfer rate M and 
corresponding period derivative P. The observed period his- 
togram however represents a sample of systems selected from 
the present-day Galactic CV population according to cer- 
tain selection effects. Thus the sample comprises systems 
with different ages at different evolutionary stages. In the 
next section we make use of the population synthesis tech- 
nique developed by Kolb (1993) and de Kool (1992) to syn- 
thesize a theoretical observable period distribution in the 
framework of a given evolutionary model for CVs. The key 
element in this procedure is to enforce the mass-radius re- 
lation as claimed by C98 on the lower main sequence. As we 
will see the resulting period distribution does indeed show 
similar spikes as the discovery probability in Fig. |l|, and no 
gap. 

3.1 A population model 

CVs are believed to form from a subset of zero-age binaries 
with a low-mass companion which undergo a common en- 
velope (CE) phase (see e.g. de Kool 1992, Politano 1996). 
The CE evolution tightens the orbit and exposes the white 
dwarf core of the giant primary. Orbital angular momentum 
losses J further shrink the orbit of the post-CE binary until 
the secondary fills its Roche lobe. The subsequent semi- 
detached evolution is governed by J and the secondary's 
response C, to mass loss, as outlined in Sect. 2. 

The method for obtaining a full theoretical period his- 
togram has been described extensively by Kolb (1993). 
Chiefly, there are three main steps involved: first, to cal- 
culate the time-dependent formation rate of nascent CVs, 
e.g. as a function of initial white dwarf and donor mass, 
for given formation rate and distributions of primary mass, 
meiss ratio and orbital separation of zero-age binaries. Sec- 
ond, to evolve this nascent population into its present-day 
configuration. The resulting multi-dimensional distribution 
function defines the intrinsic present-day Galactic CV pop- 
ulation. Third, selection effects have to be defined which 
determine how the observable sample is drawn from the in- 
trinsic population. 

In C98's attempt to explain the CV period gap the na- 
ture of the angular momentum loss J is irrelevant as long as 
it is continuous and smooth along the evolution. Obviously, 
for a population synthesis we have to specify J. We adopt 
the simplest choice for a continuous J and assume that it is 
given by the emission of gravitational waves only. For this 
case the time-dependent CV fomation rate has already been 
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Figure 2. Bipolytrope calibration functions (solid: entropy jump 
h; dashed: core index ni), used to mimic the mass-radius relation 
claimed by C98. 



calculated by Kolb & de Kool (1993), using a Monte Carlo 
technique (see also de Kool 1992). In their simulations the 
primary and secondary of zero-age binaries form indepen- 
dently from a Miller & Scalo type IMF (Miller & Scalo 1979), 
the distribution of initial orbital separations a is fiat in log a, 
and the common envelope efficiency is 1. 

To calculate the secular evolution of CVs we use a sim- 
plified stellar structure code that describes the donor's ra- 
diative core (if present) as a polytrope with index ni, the 
donor's convective envelope as a polytrope with index 1.5. 
The bipolytrope model has two free parameters, the core 
index n\ and a discontinuity h of the specific entropy at the 
stellar surface (see Kolb & Ritter 1992). We calibrate these 
parameters as a function of donor mass (Fig. ^) in such a 
way that the corresponding (bi)polytrope models in thermal 
equilibrium reproduce the C98 mass-radius relation to bet- 
ter than 1% in the range 0.12 < M2/M0 < 0.65. For larger 
masses we used the same calibration as in Kolb & Ritter 
(1992), for smaller masses a calibration which reproduces 
the CV minimum period at ~ 80 mins. To mimic the effect 
of nova outbursts we invoke an isotropic wind from the white 
dwarf, which removes the transferred material with the pri- 
mary's specific orbital angular momentum from the system. 
We note that this stellar evolution code self-consistently 
takes into account the donor's deviation from thermal equi- 
librium as a result of mass transfer, i.e. does not force C, 
to be equal to the equilibrium index along an evolutionary 
sequence. As J is rather weak the actual deviations are of 
course small. 



3.2 Results 

In Fig. ^ we plot the period histogram n(log P) of a mvis- 
limited sample taken from the resulting intrisic model popu- 
lation. To do so we neglected the drop of the CV space den- 
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Figure 3. Period histogram (liatched; scale on tlie left) and cu- 
mulative distribution (solid curve; scale on the right) for observed 
CVs with P < 5 hr. Data taken from Ritter & Kolb (1998). 
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pronounced peaks at P ~ 2 hr and 3.4 hr. The population 
density n between these spikes is not lower than outside of 
them, i.e. there is no gap. The drop of n for P ^ 3.5 hr is 
from a drop in the formation rate at larger donor mass, and 
is also seen in the period distribution obtained by Kolb & 
de Kool (1993) with the more conventional underlying main- 
sequence mass-radius exponent ~ 0.8 ~ constant. The third 
spike at the minimum period ~ 1.3 hr occurs for the same 
reason as the other two spikes, i.e. because ( — +1/3. It 
indicates the transition from the lower main sequence to the 
brown dwarf regime where —1/3. It is well-known that 
the observed period histogram (Fig |^ has no such spike 
at the lower period cut-off. This suggests that additional 
selection effects prevent us from detecting these extremely 
low-mass donor CVs. 

Changing the exponents of Mi and M in s, as would 
be appropriate for samples that are not purely rrivis hm- 
ited, does not alter the qualitatitive features of the resulting 
period histogram. Hence the detailed population synthesis 
fully confirms the expectations from the functional of the 
discovery probability (|lo|): the mass-radius relation claimed 
by C98 does generate a gap in the period distribution of 
CVs. 

From (|l^) it is easy to see that a real period gap - 
a region where p(logP) is low or zero - requires a very 
large effective (. This in turn means that we require the 
secondary's radius to change significantly for a very small 
change of its mass. This is indeed precisely what happens 
in the conventional explanation of the gap: R2 shrinks from 
a non-equilibrium value towards the main-sequence radius 
appropriate to the secondary's current mass AI2, with no 
change in the latter. Since P oc {Rl^/ M2)^^^ this radius 
change implies that the period must change if the star is 
to fill its Roche lobe, and thus transfer mass. It is hard to 
see how any revision of the equilibrium mass-radius rela- 
tion, as proposed by C98, could have the required property 
of producing a very large ^ for systems in the period gap. 

The fact that features of the function (^{Ah) describ- 
ing the lower main sequence should reappear as a mod- 
erate anomaly in the observed period histogram was al- 
ready noted by Kolb (1993). The theoretical main sequence 
has a well-established hump in ^(-^'^2) in the mass range 
0.4 ^ M2/MQ ^ 0.7, caused by the dissociation of H2 
molecules in the outermost stellar layers. This should pro- 
duce a characteristic plateau and shoulder in p(logP) at 
P ~ 5 hr. This period range is not shown in Fig. M; the data 
are still too sparse at these periods to confirm the reality of 
this feature. 



Figure 4. As Fig. |3| but for a visual magnitude— limited sample 
drawn from a model population obeying the equilibrium mass- 
radius relation claimed by C98. 

sity perpendicular to the Galactic plane and assumed that 
the visual fight comes from an optically thick, steady-state 
accretion disc. In this case the integration of the intrinsic 
space density distribution over the observable volume can 
be approximated by a multiplication with the differential 
selection factor s oc L^^ oc Mf^^M^/* (cf. Kolb 1993, Kolb 
1996). 

As expected, the period distribution in Fig. El shows two 



4 IS THE PERIOD GAP REAL? 

We have so far assumed that the period gap is genuine. 
Verbunt (1997) argues that that this may not be the case, 
essentially because not all CV subtypes show a significant 
gap. In particular, almost all (non-magnetic) novalike sys- 
tems have periods above the gap, while dwarf novae populate 
both sides of it. Verbunt concedes that the period gap is sta- 
tistically significant for the dwarf novae, but he argues that 
the gap might not be significant if one broke up the dwarf 
novae into subtypes such as the SU UMa and U Gem sys- 
tems. Almost all of the SU UMa systems are below the gap. 
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the U Gem systems above (for the period distribution of all 
CV subtypes see the catalogue of Ritter & Kolb, 1998). 

A critical re-evaluation of Verbunt's arguments has al- 
ready been given by Hellier & Naylor (1998). They convinc- 
ingly dismissed his classification of VY Scl stars with dwarf 
novae and pointed out the importance of treating magnetic 
systems separately for a comparison of dwarf novae and no- 
valike systems. Here we add to this by emphasizing that 
Verbunt's argument actually works positively in favour of 
the conventional theoretical framework of the explanation 
of the period gap. For the most basic feature of this ex- 
planation is that average mass transfer rates above the gap 
should be significantly higher than those below. Thus be- 
low the gap we expect that these rates will be too low to 
ensure the stability of the accretion discs in non-magnetic 
systems, since the disc edge will be too cool to keep hydro- 
gen ionized (cf Kolb, 1996). The prevalence of dwarf novae 
below the period gap is thus a natural consequence of the 
interrupted-braking picture and disc instability models (see 
e.g. Meyer-Hofmeister & Ritter, 1993, for a more detailed 
discussion) . 

The coexistence of dwarf novae and novalike (steady) 
systems at similar periods above the gap shows that there 
is evidently a range of mass transfer rates at a given period. 
The origin of this range is unclear, but it may result from 
long-term cycles of the mass transfer rate about the sec- 
ular mean, possibly caused by the weak irradiation of the 
secondary star by the primary (cf King et al., 1995, 1996). 
The mean mass transfer rates below the gap are so far be- 
low the critical rate for disc instability that such fluctua- 
tions evidently have amplitudes too small to produce many 
steady systems below the gap. However above the gap the 
mean transfer rate and the critical rate are quite close (cf 
Kolb, 1996), so the fluctuations can carry systems across the 
latter. Thus both outbursting (dwarf nova) and persistent 
(novalike) CVs are found above the gap. 

More generally, since the mass transfer rate —M2 and 
the mass ratio AI2/M1 are the two most important parame- 
ters characterizing the immediate appearance of a given CV 
system, it would be surprising if the explanations for vari- 
ous types of CV behaviour did not involve them. Since both 
quantities are systematically smaller below the gap than 
above, one would thus expect a preference for one side of 
the gap or the other for various subtypes, as is actually seen 
in dwarf novae. 

In any case, Verbunt's argument cannot be regarded as 
compelling, since it signally fails to provide a reason other 
than unspecified "selection effects" for the very strongly sig- 
nificant presence of a period gap when one takes the full set 
of systems in which a white dwarf accretes from a low-mass 
companion, i.e. the full CV sample: compare Fig. 0. 
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5 CONCLUSIONS 

We conclude that the conventional explanation of the CV 
period gap in terms of interrupted orbital braking is still 
the most plausible. In conjunction with other ideas, such as 
the disc instability picture, it offers a simple framework for 
understanding the period distributions of the various CV 
subtypes. This cannot be said of any proposed alternative. 
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